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Abstract. - We introduce a new family of models for growing networks. In these networks 
new edges are attached preferentially to vertices with higher number of connections, and new 
vertices are created by already existing ones, inheriting part of their parent's connections. We 
show that combination of these two features produces multifractal degree distributions, where 
degree is the number of connections of a vertex. An exact multifractal distribution is found for a 
nontrivial model of this class. The distribution tends to a power-law one, 11 (q) ~ q~"' , 7 = 
in the infinite network limit. Nevertheless, for finite networks's sizes, because of multifractality, 
attempts to interpret the distribution as a scale-free would result in an ambiguous value of the 
exponent 7. 



Networks of various kinds, such as the World Wide Web, citation network of scientific 
papers, social networks, neural networks, etc. (see are very popular objects of studies 

nowadays. There crucial role is played by the degree distribution function (DDF) 11 (q), where 
the degree q is the number of connections of a vertex (sometimes it is called connectivity). 
Data, obtained from the observations of many existing networks were interpreted as if they 
are scale-free, that is degrees of their vertices were distributed according to 11 (g) ^ |5|. 
Such networks have some specific properties, as compared, e.g., to classic random graphs |9|, 
where degree distribution follows a Poisson's law. In particular, they are extremely resilient 
to random breakdowns pO|"p^. 

To describe a scale-free growing network, the mechanism of preferential linking ^ was 
proposed. This principle is similar to the one introduced in the well known Simon model p3[ |, 
used to explain power-law distributions in various social and economic systems. In fact, all 
these models belong to the class of stochastic multiplicative processes . Here we show that 
a natural generalization of studied models of evolving networks leads to multifractal degree 
distribution with much readier properties than a power-law one. 

In general, as network is growing, two parallel processes take place: (i) New edges are 
formed between vertices. They attach preferentially to vertices of a growing network with a 
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high number of connections, (ii) New vertices appear. Several types of such a models were 
proposed [|5|,|8|,p^|-]l7t, that differ in details of how new vertices and edges are entered into the 
network. They produced scale-free networks with the 7 exponent of the DDF, either in the 
range (2, 00) [^,^, which means, that the average degree of a vertex remains finite, or 
7 < 2 ||r^,|lj], when number of edges is growing faster then of vertices. 

The questions are: whether all the scenarios of the evolution of networks, produced by 
preferential linking, lead to scale-free DDF? What other types of DDF do arise? These 
questions were considered in several recently published papers |p|,p^[- p7|p^ . However, in most 
of the models, where the idea of preferential linking was used, it was assumed, that new vertices 
appear with the same properties, independent of the state of network at this moment, i.e., 
all vertices are born equal. (Nevertheless, see the counter-example of a "breathing network" 
in Ref. One could say, that new vertices are created by some external source, whose 

properties are independent of the network's current state. Here we put forward a different 
concept: new vertices are born with random properties, which reflects the state of the network 
at the moment of birth. In this respect, one could say, that they are created by the network 



Our model may be formulated as follows: (i) New edges are introduced between existing 
vertices (see Fig.^). This happens with probability m per unit time. We assume the linear 
preferential linking: the probability, a new edge to point to a vertex i of in-degree qi is 
Pi — Qi/Q- Here Q {t) = q.^ is total degree (number of edges) of the network. We shall look 
only at the in-degrees (numbers of incoming edges) of vertices, putting aside the question of 
what vertices these edges come from. For example, one may assume that adge edge comes 
from any vertex with equal probability, (ii) New vertices appear with probability n per unit 
time (Figj^b) (in the following we assume n = 1, which can always be done by the proper 
choice of the time scale). We suppose that every new vertex has a parent — some vertex, 
randomly chosen among N (t) existing ones. Degree of a daughter vertex is assumed to be 
a random number, whose distribution depends upon parent's degree. Namely, we assume, 
that every edge, pointing to the parent with some probability c is inherited by the daughter. 
We emphasize, that here inheritance means copying of edges, — parent loses nothing, simply 
some new edges are created, pointing at its heir. If parent's degree is q, the probability for its 
heir to have the degree qi < q is {^Jc''^ (1 — c)''"'^ . 

It is the way of introducing new vertices into the network, that makes our model differ 
significantly from previously studied ones. Usually it is assumed, that new vertices have some 
preset properties, which are independent of the current state of the network. For example, 
in most cases new vertices are entered with some given initial degree qo (see Fig.|^). In our 
model new vertices are being born by network's components themselves, and with properties 
(initial degrees), reflecting the current state of the network. 



DDF may be introduced as: H [t, q) = ni/N{t)] Y^^^f 5k [q-j (t) -q]), where Sk [q - q') is 



the notation for Kronecker's (5-symbol. Averaging is over the statistical ensemble of networks, 
whose evolution is governed by the above described rules. Here N {t) is an integer, growing 
in time in a random fashion, — introduction of new vertices (as well as new links) may be 
considered as a random process (e.g., Poissonian one). Its average is iV = J dtn{t) ^ t. For 
t ^ 1, this function may be shown to obey the following master equation: 



itself. 



d_ 

dt 




OC 



^-'n(t,g') 



(1) 



Here q{t) = Q {t) /tis the average degree of network's vertices. The above equation may easily 
be derived, but, as one can see, it is a direct generalization of the master equation, introduced 
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in |T5j. The only difference is that here in the rhs, instead of 5k (q — Qq), which corresponded to 
vertices, born by some external source with definite degree, (Fig.0c), we have a term, reflecting 
DDF of a network, reached at the given moment. All terms in Eq.(^ have transparent physical 
meaning, reflecting balance between income and outcome within the number of vertices with 
degree q, tH {t, q) (note, that t variable is rather the total number of vertices in the network, 
then the time). From Eq.(|l|) one can easily obtain: q{t) = m/[{l — c)] (l + foi"^"^), where b 
is a constant of integration. Then, in the limit of large size m/q{t) may be replaced with 
1 - c. 

Coefficient c reflects the "succession right" of a given network. In principle, it can 
also be a random number, characterized by some distribution density h{c). The general 
evolution equation acquires much simpler form, if to pass to a new representation of the 
DDF, which is a bit modified Z-transform: H (i, g) -> $ (t, y) = X;^o ^ (t, q) (1 - = 

\\/N{t)\ Yli'j:^! (1 ~ vY' /■ Then the evolution equation becomes: 



{t,y)-{l-c)y{l- y) ~ j^^^ ^^^^ = ^ " 



It is to be supplied with the initial condition: <I>(to,y) = (j/), 1 (recall, that the 

equation is valid for t ^ 1). After rescaling of the size variable t — > t/to the initial condition 
becomes $ (1, y) = $o (y)- 

General solution of Eq.(|l|) may be found in the continuous approximation, when it is 
assumed, that 11 (t, q) is a slowly varying function of q. Then this equation, taking into 
account the above introduced randomness of c, takes the form: 

{t, q) + {l- c) it, <l)- I'^h (c) H (i, ^) = . (3) 

It can easily be solved after transition to Mellin's representation with respect to q: S {t, ^) = 
dq n {t, q) q^~^. Note, that S (i, n + 1) = Af„ (t) are the moments of the distribution. The 
solution is: 

s(^,a = So(0^"^«-'^ r(0 = (i-c)e-i + x(e + i) , X(0= I" dch{c)S~\ (4) 







One can see, that moments of the distribution scale with network as M„ '--^ f^^"), where 
(n) is a nonlinear function of n. These distributions are usually referred to as multifrac- 



tals 1 21 -23 1, as distinct from fractal ones, where r (n) depends on n linearly, t (n) = (n — 1) D, 
D is called fractal's dimensionality. Such a distributions were found in many objects, ranging 
from localized excitations in disordered solids to the distribution of matter in the Universe. 
Mutifractals may be thought of as a statistical mixture of fractals with different dimension- 
alities. Indeed, if to introduce the f (a) -spectrum of dimensionalities a, then a mixture of 
a-dimensional fractals with statistical weights g (a) t~^^°'\ ascribed to a-dimensional fractal, 
yields the distribution moments M„ — J da g (a) t"^"~^''~^^"\ At large t they scale with 
network's size as Af„ t'^^'^\ The prefactor g (a) gives only the coefficients of proportional- 
ity, and is not relevant. Functions / (a) and r (n) are connected by the Legendre transform: 
T (n) -I- / (a) = an, n = df /da, a = dr/dn. 

As an example, let us consider the case of the homogeneous distribution of the inheritance 
coefficient c, h{c) — d{c)9{l~c). The motivation for such a choice is that, as we shall 
see below, in this case an exact solution is possible without continuous approximation. We 
have X (C) = l/^i ^ind r (^) = £,/2 — 1 -I- 1/ (^ -|- 1). This gives multifractal / (a)-spectrum 
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to be: /(a) = (1 — -^1/2 — g)^, where — oo < a < 1/2. Note, that it includes negative 
dimensionahties (see |24 



If we set the initial condition for Eq.(|3|) as H (1, q) — Hq (q) = S {q ~ qq) , then the solution 
is a Green's function 11 (t; q, qo), whose expression is: 




(5) 



where /i is the modified Bessel's function. If t is large enough, the above formula everywhere, 
except the tail region q « qo^/t, may be replaced by: 



n (, ,„) . -^1^ exp [2f;^:^^^^ . (6) 
At t ^ oo this expression is asymptotically equal to: 

U{t;q,qo)^^{^) (7) 

q^^/^\qj {\ntf'^ 

However, this latter result is valid only if \n{q/q^)) <^ \/lnt. At t oo this region is small 
as compared with the one of the validity of Eq.(|g), which is Inq < (1/2) Int. This means 
that, in spite of that the distribution asymptotically assumes a scale- free form, ~ q^'^ with 
7 = \/2, at finite network's sizes this is true only within rather restricted region of degrees, 
small compared with the upper cut-off qoVt- However, it is extremely important, that one 
has to compare not the sizes of the regions themselves, but their logarithms. This ratio is of 
the order of l/Vlnf. Even for the largest known network, the World Wide Web, whose size is 
~ 10^, this value is of the order of 0.2. These features make the analysis of the distribution in 
terms of scale-free functions ambiguous, the value of the exponent 7 becomes dependent on 
the choice of the degree region, from which it is extracted. The (negative) slope of a log-log 
plot of the distribution steadily grows with In q from \/2 at In g < \/lni, until the exponential 
cut-off is reached at In g ~ Int. 

There is another peculiar feature of the distribution (||), obtained in the continuum ap- 
proximation. Namely, for any > at large enough t the fraction of vertices with q > qc 
scales with t as t-3/2+\/2^ This means, that the distribution concentrates within the region 
of small g, where continuum approach is invalid. Fortunately, for the homogenous distribu- 
tion of inheritance coefficient c, the exact solution is possible. Indeed, for the homogeneous 
h{c), after application of the operator dy{y-), Eq. (H), may be reduced to a linear partial 
differential equation. This equation, after the Mellin's transformation with respect to time, 
■i/)(ry, 1/) = dtt^^^^{t,y), takes the form, 

y\l-y) dli, + y (277 - Sy) dy^, + 2^^, =-29,(y$o), (8) 

where $0 (?/) is the initial distribution in the Z-representation. For the Green's function it 
is: $0(2/) — (1~2/)*'- Let us denote Mellin's time-transform of the Green's function as 
^ ^^'^q^q'^)■ Eq. (|^) may be reduced to an inhomogeneous hypergeometric one after the 
substitution -0 (y) = y'^x {v)i where C is one of the roots of the characteristic equation: Q — 
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(1 — 2ri) ( + 2ri — 0. Here we present the result in terms of 5* (?/; q, k), which is the q-th term 
of Taylor's series of ijj (y) around the point y ^ I. After lengthy calculations, we obtain 



*(?7;g,fc) 



-^i^^if^§^02(Ci,9) , 9> fc>0; 



* (77; 0, k) = -</)! (Ci, fc) ; ^ iv; q, 0) - , (9) 

where is the root of the characteristic equation, which is positive for Re 77 < 0, C2 is the 
other root, and functions (/)i,2 may be expressed as: 

, , 4gor(C)r(2 + 0/(1-0 



r(i+C + TTc)r(2+C-T^) 
, , r(i + c + T^)r(^-c)sinK 



r(T^-i)r(T^ + i) - 



(11) 



Here 2F1 is the hypergeometric function. 

Two formulas may be derived from Eqs.(p[-|ll|), which indicate: i) that the results of the 
continuous approximation are valid for large q, except for some minor corrections, and ii) that 
concentration of the distribution at low degree values is to be interpreted as in the large size 
limit, almost all vertices have zero degree. In the large q, large t limit, the following expression 
may be written for the Green's function: 



n {t; q, qo) w uj'g (qo) 



In (aq) 



(tint) 



3/2 



exp 



2yintln (yi/q^ 



9 (qo) = 



d(t)i (0 qo) 



dC 



C=%/2-l 



where: w' = 0.0823 ... and a = 0.840 ... It differs from Eq. 



(12) 

in some details, but the main 



conclusions remain the same. The formula (^ is reproduced, apart from the appearance of the 
additional multiple Ing/ Int and from different numerical factors. Also, the dependence on the 
initial degree value go is different (one may reproduce completely continuous approximation 
result (H) from the exact solution, if to assume the limit q ^ 00, qo 00, q/qo and t fixed). 
Anyway, apart from some power of logarithm, DDF values for any g > decay as /;~3/2+\/2 
with network's size increase. For g = 0, we have: 



H(t;0,go)«l-5(9o)' 



-3/2+V2 



One can see from Eq. (|l3|) that at long times the fraction of zero-degree vertices tends to 
1. These vertices do not have incoming edges but only outgoing ones (recall, that degree, in 
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the present paper, is defined as in-degree, the number of incoming edges). They are passive 
constituents of the network, i.e., their degree remains unchanged aU the time. Although the 
fraction of active vertices (with non-zero degree) tends to zero as the network grows, their 
total number increases with time as t^~^/'^ / Iji^^^ t. One can introduce the DDF of active 
vertices. Hi {t;q,qo). It follows from Eqs. ( |l^ ) and (p^), that in the large size limit, it tends 
to the distribution. Hi {q,t), which does not depend on the initial conditions: 



Hi (t, 9) ^ (1 - 6,0) , "n.^'n"^ , - t^^^^ln {aq) exp 
i — 11 (t; U, qo) 



2Wln<ln Vt/q 



(14) 



uj = 0.174 . . .. One can obtain from Eq. ( [T^ ) the following asymptotic expression at t ^ 00: 

ni(g)--^ln(ag) . (15) 

Apart of the logarithmic factor, this is a power-law (scale-free) DDF with exponent 7 < 2. 
However, the range of validity is small, Ing <C Vlni, and the moments of the distribution are 
defined by its behavior outside "scale-free" region. More general expression (|l^ ) is valid if 
q < qoVi, at (7 « qoVi there exist an exponential cut-off, which replaces the (5- functional term 
at the point of abrupt cut-off, found in continuous approximation (Eq.(H)). 

In conclusion: natural generalization of network's evolutionary dynamics with preferential 
linking was suggested. We have found that evolving networks with preferential linking, in 
which new vertices are born by previously existing ones, and inherit edges from their parent, 
arrive at the multifractal degree distribution. For a specific model, exact large-size solution was 
found. Degree distribution tends to nearly power-law type in the infinite network limit, Il{q) 
q~'^ln.q, 1 < 7 < 2 (specifically, 7 = \/2 for the exactly solvable model). Nevertheless, for 

finite network the region q < exp (^\/hi tj where this behavior may be observed, is rather small 

as compared to one of relevance, q < qoVt. With probability close to 1 a vertex, randomly 
chosen among ones with nonzero degree (active vertex), has its degree within the power-law 
dependence region. However, distribution's moments will be defined by the broad large- 
degrees tail, containing a small fraction of vertices. Results, obtained here for the model with 
a particular (homogeneous) distribution of inheritance coefficient, may easily be generalized 
to a more generic case of any distribution. 

In most of empirical scale-free distributions of degrees the exponent 7 was found to be 
7 > 2. Also, in most models, 7 > 2, that means that the average degree remains finite as 
network grows. Note that the degree distributions with 7 < 2 were obtained analytically 



and by simulation for networks with accelerating growth |17 , where q 00 as t 00. In 
the present case, the number of edges per active site also grows with time in a power law 
manner, i.e., as 

^3/2-V2 1jj3/2^^ However, the case of multifractal distribution is essentially 
different from a scale-free one. Here the distribution can not be described, using one, or finite 
set of scaling exponents. In particular, the attempt to analyze the distribution in terms of a 
scale-free one would lead to the value of the exponent 7, dependent on the method of analysis. 
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Fig. 1 - Elementary processes, defining tlie evolution of a network, (a) — creation of new link, 
pointing at a node. Probability of this process per unit time is mq{t)/Q{t). (b) — creation of new 
node by a randomly chosen member of the network. Its probability is W (q, q') = n(J,)c'' (1 — c)''~'' . 
(c) — creation of new node by external sourse. This process has the probability n. 
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